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ABSTRACT 



This study uses a numerical model to investigate the effect 
of the air-sea fluxes of sensible heat and moisture in extra- 
tropical cyclogenesis over the open ocean. The model is a 
sectorized version of the UCLA general circulation model in- 
cluding the Ara>awa-Schubert cumulus parameterization scheme 
and the Randall version of the Deardorff planetary boundary 
layer parameterization. Idealized initial conditions are speci- 
fied in the atmosphere and ocean that are typical of open-ocean 
fall and spring conditions. 

Adiabatic and diabatic model results are compared over a 
15-day integration period. Diabatic processes, including the 
surface fluxes, cause a large reduction in low-level static 
stability during the initial cyclone growth period. This re- 
duction in static stability, as well as the latent heat release, 
leads to the rapid growth of wave numbers 12 and 18 in the dia- 
batic model experiment, while only wave number 6 is present in 
the adiabatic experiment. The growth of the cyclones is much 
more rapid in the diabatic experiment. However, cyclones in 
the adiabatic experiment attain similar maximum intensities as 
in the diabatic experiment, and undergo an analogous decay period. 
After development of the initial cyclone, a variety of secondary 
low developments occur in the diabatic experiments that are 
similar to those observed over the open ocean. 
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The role of the air-sea fluxes is studied by selectively 
removing the surface fluxes at various stages of the cyclone 
development. Removal of the air-sea fluxes produces signifi- 
cant changes in the evolution of the diabatic model cyclones. 
Removal of the surface sensible heat flux results in an in- 
tensified low-level temperature gradient in the developing 
cyclone. The enhanced low-level temperature advection in- 
creases large-scale lifting and therefore latent heat release, 
resulting in an intensification of the cyclone. Removal of 
the surface moisture flux prevents the development of second- 
ary lows which might otherwise eventually replace the primary 
cyclone, as occurs in the complete diabatic model. These 
model results are not verifiable with observations, however, 
they show a mechanism by which some disturbances over the 
ocean could develop while other seemingly similar disturbances 
would not. 
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INTRODUCTION 



I . 

Extratropical cyclones are responsible for the meridional 
transport of heat and momentum as well as for most signifi- 
cant weather in the mid-latitudes. Simpson (1969) estimates 
that 75% to 85% of the mid-latitude transfer of heat, momentum 
and moisture through the ocean surface occurs during the pas- 
sage of extratropical cyclones. Understanding the general 
circulation of the atmosphere and oceans, therefore, requires 
an accurate knowledge of these large-scale eddies and how 
they utilize the large amounts of energy transferred to them 
from the world oceans. These storms also cause considerable 
destruction (Gyakum, 1980; Rasmussen, 1979) and disruption of 
commerce. Current numerical models frequently fail to pre- 
dict these systems accurately beyond one to two days, possibly 
due to the inadequate parameterization of the diabatic proc- 
esses (Sanders and Gyakum, 1980) . The success of medium and 
extended range forecasting may depend on the proper under- 
standing and modelling of the role of the air-sea fluxes in 
the development of the extratropical cyclones. 

Extensive circumstantial evidence indicates that the air- 
sea fluxes play a major role in ocean cyclogenesis. Petterssen 
et al . , (1962) described the evolution of the cyclone based on 

a composite of 51 cases of western North Atlantic Ocean storms. 
They found that initial storm development over the ocean ap- 
pears to be due mainly to low-level thermal advection. This 
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is in contrast to development over land where upper-level 
vorticity advection is the dominant development mechanism. 

Nitta and Yamamoto (1974) studied open ocean cyclogenesis 
over the East China Sea, They also found cases of cyclogene- 
sis occurring with zonal flow aloft, which implies that low- 
level thermal advection is the dominant development mechanism 
over the ocean. Other case studies (e.g., Winston, 1955; 

Pyke, 1965; Sanders and Gyakum, 1980) of explosive cyclogene- 
sis over the oceans suggest that air-sea exchanges are respon- 
sible for the observed rapid rates of development. More 
indirect evidence has been obtained recently from studies of 
cyclogenesis over the southern oceans. In a review of air- 
sea interaction in the southern hemisphere. Baker (1979) cites 
studies that show the highest frequency of cyclone development 
occurs over the oceanic polar front. Carleton (1981) obtained 
the same result using satellite data to locate the initial 
positions of developing storms in the southern hemisphere. 

These studies are only able to postulate that air-sea exchanges 
are important in over-ocean cyclogenesis. Sanders and Gyakum 
(1980) support this observation and also find a positive corre- 
lation between the intensity of the cyclogenesis and the 
strength of the ocean temperature gradient. 

Conflicting evidence regarding the importance of air-sea 
fluxes has been reported from several numerical model investi- 
gations and one observational study. Spar and Atlas (1975) 
investigated the response to sea-surface temperatures using 
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the Goddard Institute for Space Studies general circulation 
model. They found the model to be insensitive to the speci- 
fication of sea-surface temperatures for both short- and 
extended-range forecasts. Danard and Ellenton (1979) studied 
the role of the air-sea fluxes based on real data predictions 
with an eight-level primitive equation model. They also found 
that the air-sea fluxes are unimportant to modelled cyclo- 
genesis along the east coast of North America. Wei (1979) 
presents the only contrary observational study. She inves- 
tigated the energetics of an East China Sea storm using the 
Air Mass Transformation Experiment (AMTEX) data. She found 
that surface sensible heat fluxes and latent heat release con- 
tribute very little to the total available potential energy 
of the system, and concludes that the diabatic processes are 
unimportant in ocean cyclogenesis. 

The mechanisms by which air-sea fluxes could enhance ocean 
cyclogenesis are also in question. Petterssen , (1962) 

suggest the modification of low-level thermal advection as a 
mechanism. Pyke (1965) and Simpson (1969) propose an indirect 
enhancement by increasing cumulus convection which, in turn, 
increases the vorticity of the cyclone. Baker (1979) suggests 
that ocean temperature gradients induce thermally direct mean 
meridional circulations, enhancing the jet stream in their 
location and thereby determining the region of cyclogenesis . 
Recent studies of intermediate-scale cyclones developing over 
the ocean reveal similar controversies. Rasmussen (1979) 
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analyzes case studies of polar lows over the eastern North 
Atlantic Ocean. His studies indicate that these lows are 
convectively driven, implying a role for the air-sea fluxes 
similar to the one of Pyke (1965) and Simpson (1969). Reed 
(1979) and Mullen (1979) presented case studies of similar 
polar lows in the eastern North Pacific Ocean. They found 
the lows are driven by baroclinic instability processes with 
the air-sea fluxes contributing mainly to the reduction of 
low-level static stability and thereby enhancing development. 
Mullen (1979) and Reed (1979) encourage model studies to 
resolve these controversies. 

The fundamental hypothesis of this study is that the 
surface fluxes of moisture and sensible heat significantly 
modify the planetary boundary layer and the lower atmosphere 
on a tim.e scale of one to four days. These effects change 
the low-level baroclinicity and static stability of the atmos- 
phere as well as the moisture availability for latent heat 
processes. These processes may account for the observed 
differences between land and ocean cyclogenesis. 

This study should be viewed as the first in a series of 
numerical experiments at the Naval Postgraduate School that 
are designed to investigate the role of air-sea interaction 
in the development of extratropical cyclones. Its purpose 
is to examine open-ocean cyclogenesis in an idealized atmos- 
phere that is in quasi-equilibrium with the underlying ocean- 
surface temperatures. The study does not include a land-sea 
boundary and is not intended to treat the more energetic 
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cyclogenesis which occurs near the coast, presumably as a 
result of extreme air-ocean imbalances. The specific objec- 
tives of this study are to: 

1. determine the contribution of the surface fluxes to 
the growth of an unstable baroclinic wave during the 
early linear growth phase; 

2. show the mechanisms through which these fluxes 
modify the growth and decay of a baroclinic wave; 
and 

3. investigate the role of these fluxes in the develop- 
ment of secondary small-scale cyclones. 

A basic assumption of this research is that the scarcity 
of conventional observations over the ocean, and the complex 
interactions of the physical processes involved, require the 
use of a numerical model. In particular, the model should 
include a sophisticated planetary boundary layer parameteri- 
zation and analytic initial conditions to avoid uncertainties 
due to incomplete or incorrect observations over the oceans. 
Previous observational studies (see Pyke, 1965; Laevastu, 
1965; Simpson, 1969) were severely hampered by a lack of both 
conventional observations and measurements of boundary layer 
fluxes required to investigate these interactions. Previous 
numerical studies (Spar and Atlas, 1975; Gambo, 1976; Danard 
and Ellenton, 1979) used either conventional observations or 
numerical models in which the parameterizations of the dia- 
betic processes were greatly simplified. This experiment 
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utilizes a version of the UCLA general circulation model that 
is being developed by the Naval Environmental Prediction Re- 
search Facility (NEPRF) for use as the Navy Operational Global 
Atmospheric Prediction System (NOGAPS) . The NOGAPS model has 
one of the most complete parameterizations of the interactions 
between the planetary boundary layer and the free atmosphere. 
This model provides a state-of-the-art tool for investigating 
the role of the air-sea fluxes. However, this is only a 
model study and must be validated when adequate atmospheric 
observations become available. 

Chapter II of this paper describes the NOGAPS prediction 
model used in this research, and the atmosphere and ocean 
conditions chosen to initialize the model. The analysis of 
the experimental results is divided into three parts for ease 
of interpretation. These are: 

1. Chapter III--the adiabatic-diabatic comparison and 
the large-scale response to the surface fluxes, 

2. Chapter IV--the initial growth phase of the cyclone 
waves , and 

3. Chapter V — the development of the secondary cyclones. 
Chapter VI is a discussion of the results of this investiga- 
tion in terms of the three objectives. 
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THE ATMOSPHERIC MODEL 



II . 

The model chosen for this study is a state-of-the-art 
operational forecast model. Extensive evaluation and veri- 
fication of the model has been performed by research groups 
at the Naval Oceanography Command (Paine, 1980) and at UCLA 
(Mechoso a^. , 1978) . The model has been chosen by the 
Naval Oceanography Command to replace their current opera- 
tional atmospheric prediction model. Any model is an ideali- 
zation of the real world limited both by parameterizations 
of real world processes and by initial conditions. By using 
a numerical model in this study it is hoped that the results 
will aid in the interpretation of future observational studies 
and also provide insight into the difficulties of current 
forecast models in handling open-ocean cyclogenesis (Sanders 
and Gyakum, 1980) . 

The following sections describe the NOGAPS forecast model, 
the method of removing diabatic processes from the model and 
the method of initializing the model. 

A. FINITE DIFFERENCE SCHEME 

The adiabatic portion of the model is described in detail 
by Arakawa and Lamb (1977), and thus only a brief review is 
presented here. The NOGAPS code for this study was kindly 
provided by Dr. T. Rosmond of the Naval Environmental Pre- 
diction Research Facility (NEPRF) and contains all modifica- 
tions made to the system through July, 1980. A smoothing of 
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the meridional wind component on the latitude rings immedi- 
ately adjacent to the poles has been added to maintain model 
stability when computing on the fine resolution grid (2° lat 
by 2.5° long) . 

The prognostic variables for the model are the horizontal 
velocity, V, temperature, T, a form of the surface pressure, 

TT , and specific humidity, q. Additional prognostic variables 
associated with the planetary boundary layer (PBL) will be 
described in the next section. The model uses a sigma coor- 
dinate system defined as 



a 




( 2 - 1 ) 



where; 

Pj = 50 mb 

^ “ ^s ” 

p is pressure and p^ is the pressure at the earth's surface. 

In this study the lower boundary of the model is open ocean. 
The top of the model atmosphere, p^, is 50 mb. No strato- 
spheric physics are contained in this version (NOGAPS) of the 
model. The vertical distribution of the six model layers, 
assuming an arbitrary surface pressure of 1000 mb, is pro- 
trayed in Fig. 2-1. The vertical velocity, a, is carried at 
the layer interfaces, with V, T and q defined at the layer 
mid-levels. Simmons and Hoskins (1976) investigated the 
sensitivity of atmospheric numerical predictions to the number 
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Figure 2-1. Vertical distribution of model large-scale 
prognostic variables. Pressure values of 
sigma levels vary with surface pressure. A 
surface value of 1000 mb is assumed in this 
figure . 
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of model layers in the vertical. They compared results from 
two-, five-, eight- and 16-layer versions of the same model. 
They found the five-layer model accurately represented phase 
speeds to wave number 16 (the limit of their study) and growth 
rates to wave number 12. For waves shorter than zonal wave 
number 12, the five-layer version overestimates growth rates 
while the eight-layer version accurately predicts phase speed 
and growth rates to the limit of their study, wave number 16. 
Six layers is adequate for this study, however a greater reso- 
lution in the lower atmosphere would be preferable for future 
studies dealing with small-scale waves. 

The variables are staggered in the horizontal (Fig. 2-2) 
such that the meridional wind component, v, is carried at 
points north and south of the center point and the 'zonal wind 
component, u, is carried at points east and west of the center 
point. The u and v components are averaged to the center (T) 
point for diabatic and frictional computations. A coarse- 
model resolution of 4° lat by 5° long was used during initial 
model development to reduce integration time. This resolution 
was increased to 2° lat by 2.5° long for the experiments des- 
cribed below to properly handle growth rates and phase speeds 
of the large-scale cyclones. The finer resolution results in 
grid distances of 222 km by 196 km at 45 °n/S, which is also 
sufficient to resolve secondary systems with a wavelength of 
1500 to 2000 km. 

The experiments are performed on a 60 degree pole-to-pole 
sector having 24 grid points east-west and 91 grid points 
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Figure 2-2. Horizontal distribution of model large- 
scale prognostic variables. ua which is 
not shown is carried at T-points . 
(Arakawa scheme C) 
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north-south. The model uses cyclic east-west boundary con- 
ditions, which allow a continuous over-ocean propagation. 

The 60 degree sector was chosen to reduce computation time 
and space requirements. Adiabatic model studies by Simmons 
and Hoskins (1978) have shown simultaneous development of 
identical waves around the globe when the model was initial- 
ized by a normal-mode, wave number six perturbation as in 
these experiments. Although this may not be strictly true 
for this diabatic model study, the differences are only due 
to small time differences in solar heating and should not 
affect the qualitative aspects of the model experiments. 

The model utilizes a leap-frog time differencing scheme 
with a time step of three minutes in the fine mesh version. 
The model diabatic package is executed every 30 minutes, and 
is followed by a Matsuno time step to control the computa- 
tional mode and to assist in assimilation of the diabatic 
effects (see Haltiner and Williams, 1980). To avoid using 
an extremely short time step to maintain computational sta- 
bility near the poles, a smoothing technique is used on the 
zonal mass flux and east-west pressure gradient terms in that 
region (Arakawa and Lamb, 1977). An additional smoothing of 
the meridional component immediately adjacent to the poles 
is required for the fine-mesh version. 

B. DIABATIC PARAMETERIZATION 

One of the assumptions of this research is that a sophis- 
ticated PBL parameterization (Table 2-1) is required to 
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Table 2-1 



Summary of the Model Features 
Model Dynamics 



60° 


pole-to' 


-pole 


sector 


Ax 


= 2.5° 


long 


(213 km at 40°N/S) 


Ay 


0 

O 

(N 

II 


lat 


(240 km) 



Cyclic east-west boundary conditions 
Staggered grid (Arakawa scheme C) 
a vertical coordinate 

6 model layers in troposphere 

Model Diabetics 

Planetary boundary layer (Deardorf f/Randall) 

diagnostic determination of surface fluxes 
variable depth 

stratus cloud layer parameterization 
Arakawa-Schubert cumulus parameterization (Lord) 
Large-scale precipitation 
Moist convective adjustment 
Radiation (Katayama/Schlesinger) 

diurnal and seasonal variation 
cloud and water modifications 
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adequately simulate the direct and indirect effects of the 
surface fluxes on extratropical cyclone development. The 
PBL treatment in this model follows Deardorff (1972) , as 
implemented in the UCLA model by Randall (1976) . It allows 
for interactions between the PBL and cumulus cloud ensembles 
and/or a stratus cloud layer at each model gridpoint. Sur- 
face fluxes of momentum, heat and moisture are diagnostically 
determined using a bulk Richardson number based on the sea- 
surface temperature and the values of V, T and q from the 
dynamic portion of the model. These values are then used to 
predict a new PBL depth and the strength of the inversion 
jumps of V, T and q at the top of the PBL. The NOGAPS PBL 
differs from that described in Arakawa and Lamb (1977) in 
that it is constrained to remain within the lowest model layer, 
i.e., beneath approximately 800 mb, to avoid the unrealistic 
PBL depths which occurred in previous versions tested at NEPRF 
and UCLA. 

Cumulus parameterization follows the scheme of Arakawa 
and Schubert (1974) as introduced into the model by Lord 
(1978) . The base of the cumulus cloud ensembles is taken as 
the top of the PBL and provides communication between the PBL 
and the higher model layers through detrainment of mass and 
moisture at the cloud top level. Condensation also occurs 
at grid points where the air becomes supersaturated, and in 
a moist convective adjustment procedure that has been added 
to remove moist convective instability between layers not 
associated with the PBL. The radiational heating computation 
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which follows Katayama (1972) includes the effects of water 
vapor and also cloud distributions predicted by the PEL and 
cumulus parameter! zations and by large-scale precipitation 
processes . 

C. MODEL VARIATIONS 

The approach used in these experiments is to selectively 
withdraw the physical processes from the complete diabatic 
model. This results in much less model shock during the first 
few hours of integration than in experiments where physical 
processes are added to the adiabatic model integration. In 
the adiabatic case, the air-sea differences and surface winds 
can become abnormally large. This can result in excessive 
surface heat, moisture or momentum fluxes as those processes 
are added to the model. In the diabatic model, the air-sea 
fluxes are in quasi-equilibrium at the time that the physical 
processes are withdrawn from the model. Thus, there is al- 
most no shock generated and the wave cyclone is only gradually 
modified . 

The experiments are based on six different model configura- 
tions. The nomenclature to be used in describing these basic 
configurations is 

1. D — the full diabatic model as described above; 

2. A — an adiabatic model created by removing the physi- 
cal parameterization package which contains all 
diabatic effects and friction, and by adding dry 
convective adjustment to maintain stability: 
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3 . F — an adiabatic model with friction represented in 
the PEL parameterization, but not allowing any heat 
or moisture transfer to the large-scale variables; 

4. D-FS — the full diabatic model without surface sensi- 
ble heat fluxes, created by setting the surface heat 
transfer to zero within the PEL calculations; 

5. D-FQ — the full diabatic model without surface mois- 
ture fluxes, created as in (4) but for the moisture 
transfer; 

6. D-FSQ — the full diabatic model without surface fluxes 
of moisture or sensible heat, created as in (4) and 
(5) . 

It should be emphasized that this is intended to be a study 
of physical processes, not of a model parameterization scheme 
or of an initialization test of questionable or erroneous 
data. These experiments start with a system in which the 
surface fluxes are relatively small. Removing the surface 
heat flux does not remove radiative sources and sinks of heat- 
ing, nor does removing the surface moisture flux prevent con- 
densation of or the redistribution of the existing moisture 
by the various cloud parameterizations within the model. 

These experiments are similar to "spin-down" experiments 
and as such they require time before the absence of the sur- 
face transfer processes is felt. 

D. INITIALIZATION 

The model initial conditions have been chosen to study 
the air-ocean fluxes during the spring and fall regimes over 
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the open ocean. The model start time has been selected as 1 
October. This makes the northern sector a fall regime and 
the southern sector a spring (April) regime with respect to 
the solar fluxes. 

The model atmosphere initial conditions are analytically 
specified following Simmons and Hoskins (1977,1978). Zonal 
jets are centered at 35 °S and 45°N with a maximum of 52 m/s 
in the southern or spring hemisphere and 35 m/s in the northern 
or fall hemisphere (see Fig. 2-3) . The wind profiles follow 
the function u(a,(|)) = UQ(a)M(y) where: 



= 

a = 
u = 

M(y) 



latitude 

vertical coordinate 
sin (j) 

sin^(TTy^) 45°N jet 

4 1247 

' sin^(Tru‘^* ) 35°S jet 



and Ug(o) is a polynomial function approximation of the 
January mean zonal wind profile at 30°N taken from Oort and 
Rasmussen (1971). This polynomial function is Uq(o) £or the 
fall profile but is multiplied by a factor of 1.5 at all 
levels to create Uq(ct) for the spring profile. 

The jet locations and intensities are intended only to 
represent conditions which could exist during those seasons. 
April jet speeds frequently exceed 80 m/s, the maximum at- 
tained in these experiments, and jet latitudes are highly 
variable (Palmen and Newton, 1969) . The mean latitudes are 
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Figure 2-3. Initial zonal wind (solid) and temperature (dash) 
profiles. Contours 10 m/s (dot-dash 5 m/s) and 



representative of spring and fall conditions as given in Oort 
and Rasmussen (1971). In addition, these jets guarantee the 
growth of baroclinically unstable waves, as found by Simmons 
and Hoskins (1977,1978). 

Temperature fields are derived using the nonlinear balance 
routine developed by Mr. E. Barker of NEPRF. The mean Janu- 
ary geopotential heights at 30 °N from Oort and Rasmussen 
(1971) are used as the reference profile to solve for the 
geopotential field from the above analytic wind field using 
the nonlinear balance equation. The initial temperature field 
is then determined from the geopotential field using the hydo- 
static equation. The reference geopotential profile is speci- 
fied at 30°S (spring regime) in the balance routine to match 
the temperature field as closely as possible to the clima- 
tological sea-surface temperature gradient. 

In order to avoid large initial surface fluxes, emphasis 
. has been placed on the specification of low-level temperature 
and moisture profiles. In these experiments, the ocean sur- 
face temperatures are independent of longitude and fixed in 
time. Ocean surface temperatures are derived from the model 
surface air temperatures, so that initially the lowest model 
layer is near neutral stability. Some adjustments are made 
to achieve both a realistic sea-surface temperature gradient 
and a field that is consistent with the simple atmospheric 
structure represented by the specified jet intensity/location 
in each hemisphere. In Fig. 2-4 the model sea-surface tempera- 
tures are compared with smoothed values at approximately 160°W 
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Figure 2-4. Model (solid) and climatological (dash) sea-surface 

temperatures. Climatology adapted from Robinson (1976) 
at approximately 160“W for April (southern hemisphere) 
and October (northern hemisphere) . 



taken from Robinson (1976) . Climatological values for the 
southern hemisphere (spring) are April northern hemisphere 
values. The lower equatorial temperatures are desirable to 
reduce the influence of the convectively-driven Hadley circu- 
lation. The specification of open-ocean conditions (no sea 
ice) is responsible for the temperature differences near the 
poles . 

Realistic moisture profiles for the marine atmosphere are 
more difficult to specify. The intent in this experiment is 
to represent open-ocean conditions, i.e., mid-Pacific or 
eastern Atlantic. Therefore, it is natural to assume much 
higher relative humidity values in the lower model levels than 
the global mean values. Initial moisture fields for the 
model were specified with a 90% relative humidity at 900 mb, 
and a linear decrease to 40% at the highest model level (70 
mb) . Fig. 2-5 is a comparison of model relative and specific 
humidities with zonal mean values from the global climatology 
by Oort and Rasmussen (1971) for the three lowest model levels. 
Model relative humidities compare favorably with climatology 
when considered as open-ocean versus global means. Mullen 
(1979) found low-level relative humidities of 80% to 85% 
over the eastern North Pacific Ocean. Specific humidities 
are significantly different near the equator and the poles 
due to differences between mean atmosphere and model tempera- 
tures in those regions. 

The low-level winds are initially less than 3 m/s. This 
helps to prevent an initial shock or imbalance due to large 
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Figure 2-5. Model initial specific humidity at sigma levels 6, 5 and 

4 (approximately 900, 700 and 500 mb) compared with clima- 
tological northern hemisphere zonal mean values from Oort 
and Rasmussen (1971) for April (southern hemisphere) and 
October (northern hemisphere) . 



air-sea fluxes. Consequently, the PBL is set to a very 
shallow depth (10 mb) with no temperature, moisture or 
velocity jumps at the top. Rather than attempt to specify 
these values, it is intended that both the low-level wind 
field and the PBL will evolve to realistic values as the 
initial perturbation develops. 

The conditions discussed above are used for all experi- 
ments with unperturbed initial conditions (hereafter referred 
to as Z for zonal) , and as the basis for creating two differ- 
ent perturbed initial conditions. The first perturbed case 
(S) is created following Simmons and Hoskins (1977) . A 
small amplitude sinusoidal perturbation is added to the 
meridional wind component at all levels prior to balancing. 
This perturbation is in zonal wave number six with maximum 
values of .5 m/s at 45°N/S, and equal to zero at the equator 
and the poles. The associated surface pressure perturbation 
is approximately 1 mb. These initial conditions are used for 
the extratropical cyclone evolution experiments . The second 
perturbed case (R) is created by inserting similar small 
amplitude perturbations into the zonal initial conditions 
(Z) after the fields are balanced. In this case, perturba- 
tions with a maximum value of .1 m/s in wave numbers 6, 12, 
18, 24, 30 and 36 are added to the meridional wind component. 
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III. LARGE SCALE RESPONSE 



The experimental results have been divided into three 
parts. This chapter describes the large-scale, baroclinic 
wave growth in the diabatic model atmosphere, and how this 
compares to the adiabatic model experiment. The adiabatic 
experiment is also verified with a similar model experiment 
by Simmons and Hoskins (1978) . This chapter introduces the 
surface flux experiments and describes their results in terms 
of large-scale cyclone development. The next chapter con- 
siders the effects of the individual diabatic processes dur- 
ing the linear growth phase (first 72 h) of the baroclinically 
unstable waves. In the following chapter, the development 
of secondary cyclone waves in the diabatic model experiment 
is discussed. In particular, the effect of the surface fluxes 
of sensible heat and moisture on the development of the sec- 
ondary cyclone waves is examined based on these model 
experiments . 

The experiments are of two basic types — those starting at 
day 0, and those initialized at day 4 using fields from the 
full diabatic model integration (Table 3-1) . The objective 
of the first group is to establish the magnitude of the model 
atmosphere response to the diabatic parameterizations . The 
second group of experiments is performed to study the influ- 
ence of the individual diabatic processes. These experiments 
employ the technique (described in Section II. C) of selectively 
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Table 3-1 



Experiments for the Large-scale Study 



* * 



Exp. 


Initial 


Integration 


Model 




Fields 


Period 


Version 




Experiments Initialized from Day 


0 


D 


S (wave no. 6) 


15 days 


complete diabatic 


A 


S 


15 days 


adiabatic 


DZ 


Z-zonal 


3 days 


complete diabatic 


DR 


R-multiwave 


5 days 


complete diabatic 




Experiments Initialized from Day 


4 


F 


S4 (fields from 


6 days 


adiabatic with 




D at day 4) 


model friction 


D-FS 


S4 


6 days 


diabatic w/o surface 
sensible heat flux 


D-FQ 


S4 


6 days 


diabatic w/o surface 
moisture flux 


D-FSQ 


S4 


6 days 


diabatic w/o surface 



moisture or sensible 
heat fluxes 



* 

See Sections II. C and II. D for a more complete discussion 
of model variations and initial conditions. 



38 



removing the model parameterizations to investigate the 
contributions of the surface fluxes of sensible heat and 
moisture . 

A. DIABATIC -ADIABATIC COMPARISON 

The adiabatic model experiment (A) compares very closely 
with the results of a similar study by Simmons and Hoskins 
(1978) . Surface pressure maps for experiment A at day 6 are 
displayed in Figs. 3-1 a-b . The northern hemisphere case 
(3-la) , which employs initial conditions adapted from Simmons 
and Hoskins (1978) , is essentially the same as in their study 
(see their Fig. 1) . The surface pressure field exhibits a 
pure wave number six pattern in each hemisphere. The growth 
of the cyclone-anticyclone pair is much greater in the spring 
(southern) hemisphere, as would be expected with the more in- 
tense initial jet in that hemisphere. Southern hemisphere 
maps in Fig. 3-lb and below are presented in inverted form 
(south at top, east remains toward right) for ease of comparison 
with the northern hemisphere results. The positions of the 
cyclone centers also reflect the more equatorward initial 
latitude of the upper-level jet maximum in the spring hemi- 
sphere. The vertical structure of the growing cyclone (Fig. 

3-3) is that of a typical growing baroclinic wave as seen in 
observational studies (Palmen and Newton, 1969) . This struc- 
ture also follows closely that obtained by Hoskins (1978) 
in a theoretical study of baroclinic waves. In agreement 
with Simmons and Hoskins (1978), the wave cyclone in A reaches 
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Figure 3-1. Surface pressure for A after six days of 

integration for a. northern hemisphere and 
b. southern hemisphere. Contour interval 
4 mb. Southern hemisphere inverted with pole 
(south) at top (east remains toward right) . 




Figure 3-2. As in Fig. 3-1 except for diabatic model 
(D) . 
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a maximum intensity at about day 9 and then decays slowly 
through day 15 of the model experiment. Simmons and Hoskins 
(1978) demonstrated that this decay was brought about by 
barotropic processes which became dominant during the later 
part of the forecast period. 

The full diabatic model experiment, D, is used for adia- 
batic-diabatic comparisons and as the control case for the 
surface flux model experiments. The inclusion of diabatic 
processes in this experiment results in much more rapid de- 
velopment of the primary low by day 6 (see Figs. 3-2 a-b 
and 3-1 a-b). In addition, shorter wavelength disturbances 
grow in the diabatic experiment as found by Haltiner (1967) . 
Although shorter wavelengths are not apparent at day 6 in 
the spring case (Fig. 3-2b) , they are present at earlier 
(see Fig. 3-6b below) and later times (see Chapter V) in 
the diabatic model run. Secondary cyclone waves are more 
prevalent in the fall hemisphere, but they are an important 
feature of both hemispheres in the diabatic model experiments. 

Experiments A and D are initialized with a wave number 
six perturbation. Following Simmons and Hoskins (1978) and 
other studies of baroclinic wave growth, it was assumed that 
wave number six would be the most unstable for these experi- 
ments. Instead, the initial period of one to three days is 
dominated by wave number 18 growth. As mentioned above, a 
wave number 6 disturbance eventually evolves. To estimate 
the sensitivity to the use of wave number 6 in the initial 
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conditions, a second diabatic experiment (DR) was performed 
using the multi-wave initial conditions described in Section 
II. D. Wave number 18 again dominates the early growth period, 
but a distinct wave number 6 pattern evolves by day 5 (Figs. 
3-4 a-b) . Haltiner (1967) conjectured that the evolution of 
large-scale waves in the atmosphere was more the result of a 
change in scale of shorter wavelengths, than the growth from 
a small amplitude perturbation of large-scale dimensions. 

This is apparently true in D and DR. This aspect of the dia- 
batic model experiments and the growth of the smaller scale 
waves is discussed in detail in the following chapters. 

The primary trough in D exhibits a larger phase speed 
than that of the primary trough in A during the period of 
strongest development (Figs. 3-5 a-b) . This speed differ- 
ence is approximately 5° long per day (450 km at 35°S) at 
its maximum; however, the phase speeds are approximately 
the same after day 5. 

The release of latent heat narrows and strengthens by 
several fold the zone of upward motion. This increases the 
conversions of eddy available potential energy to eddy kinetic 
energy, and enhances the growth of the baroclinic wave (Palmen 
and Newton, 1969; Haltiner, 1971) . This more rapid growth 
is evident in D. A result not anticipated from previous 
work is that the cyclones in A eventually (about day 9) be- 
come nearly as intense, 980 mb/980 mb (northem/feouthern 
hemisphere) central pressure versus 982 mb/976 mb, as in D. 
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Figure 3-4. Surface pressure for DR after five days 

of integration for a. northern hemisphere 
and b. southern hemisphere. Contour 
interval 4 mb. 
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Figure 3-5. Position of streamline troughs at sigma 

level 3 (300 mb) in the southern hemisphere 
for a. A and b. D. 
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Figure 3-6. As in Fig. 3-2 except after four days of 
integration. 
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Beyond day 8, the primary lows begin a gradual decay in 
intensity as in A. The barotropio decay process appears to 
control the diabatic experiment as in the adiabatic experi 
ment. Thus, although the minimum surface pressure appears 
to depend on the large-scale baroclinic and barotropic proc- 
esses, the diabatic processes greatly affect the timing and 

rate of development of the primary lows. 

in summary, the adiabatic model experiment is consistent 
with previous studies, notably Simmons and Hoskins (1978) . 

The primary cyclones attain approximately the same maximum 
intensity in the diabatic as in the adiabatic experiments 
and undergo a similar decay period. The growth of the pri- 
mary cyclone is much more rapid in the diabatic experiment, 
however, and involves interactions with wave scales not present 
in the adiabatic experiment. Also, in the adiabatic experi- 
ment the hemispheric differences are primarily m the inten 
sity of the primary cyclone. The diabatic experiment demon- 
states a stronger tendency toward shorter wavelengths in the 
fall hemisphere as well. These results agree closely with 
observational studies. Sanders and Gyakum (1980) found a 
greater frequency of rapidly deepening storms over the oceans 
where diabatic processes are more dominant. Also, there is 
a greater tendency for shorter wavelength cyclones to deve p 
over the mid- and high-latitude oceans than over land (Reed, 
1979; Nitta and Yamamoto, 1974). 
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B. SURFACE FLUX MODIFICATION TO THE PRIMARY LOW DEVELOPMENT 



Four model versions are initialized with fields from day 
4 of D (Figs. 3-6 a-b) to study the contributions of the 
diabatic processes (Table 3-1) . The broad areas of low pres- 
sure on the right side of each figure evolved into intense 
cyclones (see Figs. 3-2 a-b) during the next two days of 
integration in all model versions. Removal of the surface 
fluxes produced relatively minor changes in the intensity of 
the primary cyclones and anticyclones (Table 3-2) and very 
slight changes in phase speed (not shown) by day 6, which is 
two days after removal of the fluxes. Removal of the surface 
sensible heat flux (D-FS) leads to an additional 1-3 mb deep- 
ening of the primary cyclone. The largest effect occurs in 
the northern or fall hemisphere with the warmer SST and weaker 
upper-level jet. The removal of the surface moisture flux 
(D-FQ) leads to 4-5 mb less deepening of the primary cyclone 
in two days. No further deepening occurred beyond day 6 in 
D-FQ and D-FSQ . Thus, there are differences of 9-10 mb in 
central pressure by day 8, four days after initialization. 
These minor changes are expected from previous observational 
studies (see Palmen and Newton, 1969) and also from numerical 
model studies (Haltiner, 1971; Spar and Atlas, 1975). 

The surface pressure fields at day 6 of the no surface 
flux experiments (Figs. 3-7 a-c) are very similar to those 
of the diabatic experiment (Fig. 3 -2a) . However, large sur- 
face fluxes of moisture (Fig. 3-8) which are present in D 
and D-FS, but are not present in D-FQ and D-FSQ, lead to 
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Table 3-2 



Central pressures (mb) for the primary cyclone/ 





anticyclone by 


experiment for 


days 4 , 6 


and 8 . 


Model 

Exp. 


Hemisphere 


Time 
Day 4 


Day 6 


Day 8 


D 


Fall 


1012/1020 


988/1028 


983/1024 




Spring 


1002/1020 


975/1022 


981/1024 


D-FS 


Fall 


initialized 
with fields 
from D above 


987/1027 


980/1022 




Spring 


II 


974/1024 


980/1028 


D-FQ 


Fall 


initialized 
with fields 
from D above 


992/1024 


992/1024 




Spring 


II 


979/1022 


985/1024 


D-FSQ 


Fall 


initialized 
with fields 
from D above 


992/1024 


993/1024 




Spring 


II 


978/1024 


985/1024 


F 


Fall 


initialized 
with fields 
from D above 


997/1029 


992/1028 




Spring 


II 


997/1026 


981/1025 
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Figure 3-7. Northern hemisphere surface pressure at 

day 6 for a. D-FS, b. D-FQ and c. D-FSQ. 
Contour interval 4 mb. 
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Figure 3-8. Northern hemisphere surface moisture flux 
rate at day 6 for D. Contour interval 
4 mm/day . 
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different cyclone evolutions in these experiments. The 

largest moisture fluxes, which are equivalent to approxi- 

2 

mately 1200 cal/cm /day of latent heat, are located in the 
southwest quadrant of the primary cyclone in D and in D-FS 
(similar to D but not shown) . These values are consistent 
with previous observational studies (e.g., Petterssen et al . , 
1962 and others) . After two days, the loss of the moisture 
flux has led to large changes in the distribution of low- 
level moisture (Figs. 3-9 a-d) . Very dry air is advected 
into the primary cyclone in D-FQ and D-FSQ (Figs. 3-9 c-d) 
as a result of the lack of the surface moisture flux. This 
results in a decrease in latent heat release in the primary 
cyclone and a large reduction in its growth after day 6 
(Table 3-2) . 

It is generally accepted that the temperature gradient 
associated with the cold front of a developing cyclone is 
greater than that of the warm front. For the low-level 
temperature fields of the experiments performed here (Figs. 
3-10 a-e) this is only true in the adiabatic experiment, F, 
and possibly the no surface flux experiment, D-FSQ. Both the 
surface sensible heat flux and the surface moisture flux re- 
duce the strengthening of the low-level temperature gradient 
by thermal advection. D-FS and D-FQ have a larger tempera- 
ture gradient in the region of the primary cyclone than in D. 
This gradient is even larger in D-FSQ, and approaches that 
of the adiabatic experiment, F. Although this does not 
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Figure 3-9. Northern hemisphere specific humidity fields 
at sigma level 6 (900 mb) at day 6 for 
a. D, b. D-FS, c. D-FQ and d. D-FSQ. 
Contour interval 2 g/kg. 
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Figure 3-10. As in Fig. 3-9 except for temperature 
fields. Contour interval 5°C. Also 
shown (next page) e. experiment F. 
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apparently affect the evolution of the primary cyclone on a 
two-to-three day time scale, it could have significant 
effects on the development of secondary or subsequent 
cyclones . 

Diabatic processes play a greater role in the growth of 
the primary cyclone in the fall hemisphere with the weaker 
jet and higher SST. Over the open ocean, the surface fluxes 
of moisture and sensible heat only weakly modify the develop- 
ment of the primary cyclones on a time scale of less than two 
days, as found in the case studies of Mullen (1979) . The 
modifications of the low-level temperature and moisture dis- 
tributions are significant, however. In particular, the sur- 
face fluxes tend to weaken the temperature and moisture 
gradients associated with strong advection, especially in 
the western sector of the cyclone where the surface fluxes 
are largest (Petterssen ^ , 1962). This results in more 

intense gradients in the eastern sector than in the western 
sector of the primary cyclone in the model marine atmosphere. 
This is in contrast to that found in the adiabatic experi- 
ments or in typical over-land examples of cyclogenesis. Be- 
yond two days, these modifications affect the development of 
the primary cyclone. The response of secondary or subsequent 
cyclones may also be affected, as presented in Chapter V. 



56 



IV. INITIAL DEVELOPMENT 



A major feature of the diabatic experiment is the evolu- 
tion of multiple lows and secondary lows (see Fig. 3-2a) . A 
cyclone scale initially much shorter than wave number six 
expected from the adiabatic experiment predominates during 
the first three to four days of the diabatic experiment. 

Beyond day 4, as seen in Figs. 3-6 and 3-2, these multiple 
lows gradually evolve into an intense primary low with asso- 
ciated secondary wave cyclones. The evolution of these second- 
ary cyclones is addressed in Chapter V. This chapter describes 
the initially linear growth phase and the reasons why the 
smaller-scale cyclones occur. 

Table 4-1 lists the model experiments viewed in this 
chapter. The initial conditions and model configurations 
used for these experiments are presented in Chapter II. 



Table 4-1 

Experiments for the initial development study 
* * 



Exp. 


Initial 


Integration 


Model 




Fields 


Period 


Version 


D 


S — wave no. 6 


15 days 


complete diabatic 


A 


S 


15 days 


adiabatic 


A36 


Si. 5 (fields 
D at day 1.5) 


from 2 days 


adiabatic 


D-FS 


S 


3 days 


diabatic w/o surface 
sensible heat flux 


D-FQ 


S 


3 days 


diabatic w/o surface 
moisture flux 



* 

See sections II. C and II. D for a more complete discussion 
of model variations and initial conditions. 
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A. DIABATIC CONTRIBUTIONS TO WAVE GROWTH 

Growth rates by wavenumber at 300 and 900 mb for experi- 
ments A, A36 and D are displayed in Figs. 4-1 a-d. The 
growth rates are obtained by performing a Fourier analysis 
of the meridional wind components for two time periods. The 
exponential growth rate for each wave number is then obtained 
from 






-1 






(4-1) 



where : 

e = growth rate 

V = meridional wind component 

k = wavenumber (6,12,18,24) 

Z = model sigma level (~ 900 mb and 300 mb) 

= initial time, and 
t 2 = final time. 

The most obvious difference between experiments A and D 
is that only wave number 6 grows in A. This has already been 
described and is most clearly seen in the surface pressure 
fields displayed in Figs. 3-1 and 3-6. Experiment A36 was 
included to incorporate any changes in the model basic flow 
that may have occurred due to slight imbalances in the initial 
conditions, or due to the adjustment of the diabatic model to 
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Figure 4-1. Exponential growth rate between days 2.5 and 
3.5 for the meridional velocity perturbation 
by wavenumber (see text for description) for 
a. sigma level 3 (300 mb), 36°S, b. sigma 
level 3, 46°N, c. sigma level 6 (900 mb), 
36°S and d. sigma level 6, 46°N. Experi- 
ments D, A and A3 6 are compared. 



59 



the underlying sea-surface temperatures. Experiments A and 
a 36 use the same adiabatic version of the model — the only 
difference is that A36 is initialized using fields from day 
1.5 of the diabatic experiment while A uses the day 0 fields. 

In A36, wave number 12 grows more rapidly than wave number 6 
in the lower atmosphere (900 mb) . At higher levels, wave 
number 12 exhibits no growth in A36. Another readily evident 
feature is the north/south variation in the A3 6 experiment at 
higher wave numbers. Wave numbers 18 and 24 grow weakly at 
low levels and not at all in upper levels in the weaker jet/ 
higher SST fall hemisphere. In the spring hemisphere, these 
wave numbers exhibit a growth comparable to wave number 6 at 
both levels. The last feature to note relative to A36 is 
that wave number 6 grows more rapidly in this experiment than 
in A. 

The diabatic experiment D, exhibits a definite shift of 
the most unstable wave to shorter wavelengths. At upper 
levels, wave number 12 is the most rapidly growing wave in both 
hemispheres , while at low levels wave number 18 dominates . 

This is a distinctly different result from that of the adia- 
batic experiments. Wave number 6 also grows more rapidly in 
D than A everywhere except at low-levels in the northern or 
fall hemisphere. Figs. 4-2 a-b give the streamline troughs 
versus height and longitude for A and D. These figures pro- 
vide a clearer picture of the actual wave number/height rela- 
tionships described in the growth rate calculations, although 
variations in intensity are not indicated. 
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Figure 4-2. Location of streamline troughs at day 3 
for a. A and b. D. 
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Several model studies, e.g., Staley and Gall (1977) and 

Blumen (1979) , have found that the wavelength of the most 

unstable wave in the atmosphere is very sensitive to changes 

in low-level static stability. In addition, Rao and Ferreira 

(1979) performed a numerical study of intermediate-scale 

(less than 2000 km) waves using a moist model with a CISK- 

type parameterization of convective heating. They found that 

the use of a standard atmosphere value for static stability 
2 -2 -2 

of . 02 m s mb produced two modes of unstable waves. The 

first was a conventional, long baroclinic wave in which latent 

heating was relatively unimportant. The second was an inter- 

mediate-scale disturbance in which latent heating-induced 

vertical velocities were of the same order or larger than 

those of the baroclinic processes. When a small value of 

2 -2 -2 

static stability (.002 m s mb ) was used, only a conventional 
baroclinic wave was found, but with an intermediate-scale 
wavelength . 

The above study and a similar one by Gambo (1976) suggest 
that there are two modes of intermediate-scale unstable waves: 

1. a latent heat-driven mode with energetics similar to 
a tropical system but not dependent on low values 

of static stability; and, 

2. a conventional baroclinic mode which arises due to 
very low values of static stability. 

Clearly only the conventional, baroclinic mode can occur in 
the adiabatic experiments. While the diabatic model does not 
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have a CISK latent heat parameterization, the more complete 
parameter! zations of convective, PBL and large-scale latent 
heating processes should allow the formation of either inter- 
mediate-scale mode. In the discussion below, differences in 
wave growth among experiments A, A3 6 and D are explained in 
terms of these two unstable modes. 

Static stability is defined as 



where ; 

a = specific volume 
0 = potential temperature and 

p = pressure. 

Figs. 4-3 a-b give the variation of 39/0p as a function of 

time for experiments A, A3 6 and D for the 900-500 mb layer. 

The corresponding values of static stability range from .015 
2 -2 -2 

to .012 m s mb . Values of static stability are computed 
separately for the 900 mb-surface (D only) , and 900-700 mb and 
700-500 mb layers. The lowest values are found in the 900- 
700 mb layer (.011 to .013 m^s ^mb , but the greatest tem- 
poral variability is found in the 700-500 mb layer (.015 to 
2 -2 -2 

.022 ms mb ) . The averaging of the two model layers 
seems to describe the variability of static stability in 
the lower model atmosphere better. 



63 



Aa/AP (*K/ lOOMB) A0/AP(*K/IOOMB) 





Figure 4-3. Change in potential temperature differences 
with pressure between sigma levels 4 and 6 
(500 to 900 mb) in °C/100 mb for a. 36°S 
and b. 46°N. 
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In the region of the jets in both hemispheres, the 
zonally-averaged static stability reduces rapidly during the 
first two to three days of integration in D. The trend then 
reverses around day 3 to 3.5, and the static stability in- 
creases to a value comparable to A by day 5. Both experi- 
ments A and A36 maintain a relatively constant static stability 
throughout the period being considered, with the A3 6 case hav- 
ing a much smaller value due to the initialization at day 
1.5 from D. It is important to note that the value of static 
stability is similar in experiments A36 and D, particularly 
in the southern or spring sector. 

Low-level static stability (Fig. 4-3) and the growth rates 
of wave numbers 6 and 12 (Fig. 4-1) appear to be associated, 
especially in the adiabatic experiments. The only energy 
source for the waves in the adiabatic experiment is baro- 
clinic instability. No changes in static stability due to 
air-sea fluxes or latent heat release can occur in the adia- 
batic experiments. Thus, it appears significant that wave 
number 12 develops in A36 at low-levels, but does not develop 
in A which has a larger static stability (by a factor of 40%) . 
Also, Mullen (1979) found that small-scale cyclones in the 
north Pacific Ocean generally have very, small values of static 
stability. As noted in Gambo (1976) and Rao and Ferreira 
(1979) , the intermediate-scale baroclinically unstable mode 
is a shallow wave, and it is not present at 300 mb in A36. 

The deeper and longer (wave number 6) wave is relatively in- 
sensitive to low-level static stability changes and its growth 
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rate in A36 increases only slightly over that of wave number 
6 in A. The variations in growth rates of wave numbers 18 
and 24 at the 90 0 mb level in A3 6 may be similarly explained 
in terms of low-level static stability. However, the growth 
of these waves at 300 mb (southern case only) is nearly iden- 
tical to that of wave number 6 in A36. Because these waves 
(18 and 24) are expected to be shallow waves, their occurrence 
at 300 mb is likely due to nonlinear effects (Williams, 1965) . 

Low-level static stability in D is not significantly dif- 
ference from that of A36. Nevertheless the wavelength of 
the fastest growing wave is shifted to shorter wavelengths 
in D at all levels. Wave number 12 dominates at 300 mb and 
wave number 18 dominates at 900 mb. However, wave number 6 
grows only slightly faster in D than in A. Static stability 
is not likely to be the cause for this additional shift, and 
the only model differences are in the diabatic processes. This 
implies a latent heating process is responsible for the en- 
hanced growth of the shorter wavelength waves, in agreement 
with the second unstable mode of Gambo (1976) and Rao and 
Ferreira (1979) . Another fact supporting this reasoning is 
the much larger growth rates at 300 mb in the northern or fall 
case. In terms of baroclinic instability processes, the much 
stronger baroclinically unstable jet is located in the southern 
case. Consequently, the largest growth rates should be pre- 
sent there, which is true for the shorter waves of A36. If, 
however, latent heat release is responsible for the larger 
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growth rates, the lower vertical wind shear and higher sea- 
surface temperatures in the northern case should produce the 
largest growth. Also, as is shown below, wave number 18 com- 
mences rapid growth at the onset of D, when the low-level 
static stability is comparable to that in A. This growth 
appears to be in response to low-level precipitation or 
latent heat release occurring due to the very moist initial 
atmosphere, and also due to surface heat fluxes. 

Fig. 4-2b also depicts a more vertical structure asso- 
ciated with the shallower wave in D than in the deep waves, 
or in wave number 6 (Fig. 4-2a) of A. The east-west tilt of 
the deeper waves is more typical of a growing baroclinic wave, 
while the lack of tilt of wave number 18 is more suggestive 
of a latent heat-driven wave (Palmen and Newton, 1969) . 

To summarize this section, much shorter wavelength waves 
are present in the diabatic experiment than in the adiabatic 
experiments. Two processes seem to be responsible for the 
growth of the shorter wavelength waves. The first is the 
reduction of the low-level static stability by the diabatic 
processes. The second is the release of latent heat. The 
growth rates in the fall case with higher sea-surface tempera- 
tures and a weaker jet responded to these diabatic effects 
more strongly than in the spring case with the stronger baro- 
clinically unstable jet. These results are in close agreement 
with the case studies of Reed (1979) and Rasmussen (1979). 

Reed observed small values of low-level static stability 
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associated with the development of intermediate-scale cy- 
clones while Rasmussen observed strong latent heat release 
and a tropical structure. Both observed characteristics are 
present in the diabatic experiment and appear to contribute 
to the growth of the intermediate- sc ale waves. 

B. SURFACE FLUX MODIFICATION OF GROWTH 

The surface fluxes during the initial growth phase are 
very small in terms of normal atmospheric values (Table 4-2) . 
This is due to the specification of the initial fields and 
also because the surface wind speed associated with the small 
amplitude perturbation is still small. However, it will be 
shown that these small surface fluxes, through interaction 
of the PBL with the model diabatic processes, have a signi- 
ficant effect on the model atmosphere. 

Table 4-2 

Surface fluxes of zonal mean sensible heat 

(S) and moisture (Q) in W/m^ at day 2. 

(Q is in latent heat equivalent.) 



Latitude 


Flux 


Experiment 








D D-FS 


D-FQ 


46°N 


S 


18.3 


25.6 




Q 


48.3 42.0 


- 


36®S 


S 


15.2 


39.5 




Q 


42.5 34.4 


- 


The static 


stability 


of the surface flux 


experiments 


decreases in the same sense as in D (Fig. 4-4 


a-b) . D-FS 
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Figure 4-4. As in Fig. 4-3. S and Q represent experi- 
ments D-FS and D-FQ respectively. 
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(no surface sensible heat flux) follows D closer than does 
D-FQ (no surface moisture flux) . Comparison with experiments 
A and A36 (Fig. 4-3) indicates that static stability in ex- 
periments D-FS and D-FQ follows that of A3 6 more closely than 
that of A. Consequently, the reduction in static stability 
should increase the growth of the shorter waves in D-FS and 
D-FQ to a level at least comparable to A36. This shift, how- 
ever, should occur gradually as static stability is reduced 
and not during the initial 24 to 48 h of integration. 

Figs. 4-5 a-b give the amplitudes of the meridional wind 
components at day 2 and day 3 by wavenumber for the diabatic 
experiments. In the preceding section wave growth rates be- 
tween days 2.5 and 3.5 are derived for the adiabatic/diabatic 
comparison. Here the actual amplitudes indicate the temporal 
change of wave growth in response to static stability and 
other factors. Note also that values given in Fig. 4-5 are 
vertically averaged, whereas the values in Fig. 4-1 are at a 
single level. The values for wave number 6 are much larger 
than the other wave numbers at day 2 because the model was 
initialized with a value of approximately .5 m/s in wave num- 
ber 6. This results in a larger amplitude for this wavenumber 
even though it may have a smaller growth rate than the higher 
wave numbers. The exponential growth rates, which depend on 
the ratios of the day 3 to the day 2 amplitudes, are not 
shown here. 

The largest scale wave, wave number 6, is relatively 
insensitive to the removal of the surface fluxes. D-FS (no 
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Figure 4-5. Amplitude (m/s) of the meridional wind per- 
turbation by wavenumber for day 2 (stippled) 
and day 3 (hatched) for a, 46°N and 
b. 36°S. D, S and Q as in Fig. 4-4. 
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surface sensible heat flux) and D-FQ (no surface moisture 
flux) exhibit 10% to 15% less amplitude in wave number 6 than 
does D. There is also a smaller growth rate in the fall case 
compared to the spring case with the stronger jet. The re- 
duction in amplitude in D-FS differs from the study of dia- 
betic effects on baroclinic instability discussed in Chapter 
III. This may be due to the cooperative reduction in surface 
moisture flux in this experiment as noted in the results pre- 
sented in Table 3-2. 

Wave number 12 is the slowest growing wave during the 
first one to two days in all the diabatic experiments. Be- 
yond this time, however, it becomes the fastest growing wave 
in experiments D and D-FS. Referring to Fig. 4-4, the corre- 
lation with the decrease in low-level static stability appears 
to be very high. 

Wave number 18 shows the greatest initial growth in the 
diabatic experiments. This wavelength and wavenumber 24 are 
also the most sensitive to the surface fluxes of sensible heat 
and moisture. The amplitude of wavenumber 18 at days 2 and 
3 in D-FS is 20% to 40% less than in D. For D-FQ, the ampli- 
tudes in both hemispheres are less than 50% of those in D. 

Wave number 18 growth is more rapid in the fall or weaker jet 
hemisphere in the diabatic experiments. This behavior is 
opposite to that of wave number 6, which grows more rapidly 
in the spring hemisphere. This difference of growth could be 
either a response to lower values of vertical wind shear as 
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found by Rao and Ferreira (1979) , or a response to higher 
sea-surface temperatures. Considering changes in growth 
rate of wavenumber 18, the greatest decrease from loss of 
the surface sensible heat flux during the first two days oc- 
curs in the fall sector. However, wavenumber 18 grows at 
the same rate between days 2 and 3 in the fall case for all 
three experiments, which implies that some other process has 
become dominant. 

The vertical distribution of the meridional velocity per- 
turbation for experiments D, D-FS, D-FQ and A is given in 
Figs. 4-6 a-b. A Fourier analysis routine is used to filter 
out the mean meridional velocity and all perturbation com- 
ponents with a wavenumber higher than 30. Several tendencies 
already noted in Figs. 4-1 and 4-5 are present here. The 
largest amplitude in A is in the spring hemisphere, whereas 
it is in the fall hemisphere for the diabtaic experiments. 

The stronger spring jet controls the growth in A. The pertur- 
bation in A is much smaller than that found in the diabatic 
experiments — as much as a factor of 10 smaller than D at 900 
mb at 46°N. Removal of the surface fluxes, D-FS and D-FQ, 
makes the amplitude of the perturbation approach that of the 
adiabatic experiment. Withholding the surface moisture flux, 
D-FQ, has a larger effect than removal of the surface sensible 
heat flux. The contributions of the higher wave numbers (see 
Fig. 4-5), which grow more rapidly in the fall hemisphere, are 
responsible for the larger amplitude in that hemisphere of 
the diabtaic experiments. 
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Figure 4-6. Amplitude (m/s) at day 3 of the meridional 
wind perturbation summed over wavenumbers 
6 through 30 for a. 46°N and b. 36° S. 
Experiments D, S(D-FS), Q(D-FQ) and A are 
compared. 
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The largest amplitude in D is at the lowest model level. 
This is due to the contributions of the high wavenumber com- 
ponents which tend to be associated with shallow systems. 

In D-FQ and D-FS , the greater response of the shorter waves 
to the loss of the surface fluxes causes the greatest reduc- 
tion in growth at the lowest level. Again, D-FS and D-FQ tend 
toward the adiabatic experiment which has its largest ampli- 
tude at 300 mb. 

The temperature perturbations are much larger in the 
diabatic experiments than in A (Figs. 4-7 a-b) . The inclusion 
of the surface fluxes increases the amplitude of the pertur- 
bations. The percentage decrease in the temperature pertur- 
bation of D-FS, D-FQ and A relative to D does not reflect the 
larger decrease in wave growth that is found for the meridional 
velocity perturbation. The implication here is that the con- 
tribution to growth of the wave is not in direct response 
to the addition of low-level heat through the surface fluxes. 
These fluxes must also indirectly cause the greater wave 
growth exhibited by the meridional velocity component. 

The vertical velocity perturbations (Figs. 4-8 a-b, derived 
as for the meridional velocity perturbation) appear to be the 
link between the temperature perturbations and the wave growth 
as seen in the velocity perturbations. The vertical velocity 
perturbation in D is 10 to 20 times that of A. Haltiner 
(1971) describes how larger vertical velocities resulting 
from diabatic heating lead to a greater conversion of eddy 
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available potential energy to eddy kinetic energy, resulting 
in enhanced wave growth. A large portion of the difference 
between vertical velocities in D and A is clearly due to the 
surface fluxes of sensible heat and moisture, with surface 
moisture flux (D-FQ) having the greatest effect. The hemis- 
pheric difference is also evident here. At 46®N with a 
higher SST and weaker jet, D-FS has a vertical velocity per- 
turbation 40% to 50% less than that in D. In the spring 
hemisphere, the difference is less than about 10%. The slightly 
larger value at 400 mb in the spring case becomes slightly 
less than in D when the values are averaged over 34° to 36°S. 

In the spring hemisphere, vertical velocities associated with 
the stronger jet may be sufficient to trigger the release of 
latent heat. In the fall hemisphere with the weaker jet, the 
contribution of the larger surface sensible heat fluxes in 
generating vertical velocities may be a major factor in trigger- 
ing the release of latent heat. Hemispheric differences are 
also evident in D-FQ (no surface moisture flux) . In the fall 
case, vertical velocity perturbations range from an approxi- 
mately 85% reduction from those in D at 200 mb to a 65% reduc- 
tion at 800 mb. The corresponding spring case range is 70% 
to one of 40%. The surface sensible heat and moisture fluxes 
appear to be cooperative (as in Table 4-2) , and the lower SST 
of the spring hemisphere results in a smaller moisture flux 
or influence by the moisture flux. 

The contributions of the surface fluxes to the distribu- 
tion of the vertical velocity is not the same for both moisture 
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flux and sensible heat flux. In D-FQ, the largest perturba- 
tion is in the lowest layer, 800 mb, in both hemispheres. 

This is due to the surface sensible heat flux causing the 
greatest heating at this level. In D-FS, the maximum is at 
600 mb, due to heating caused by the release of latent heat 
at higher levels. Again, a cooperation between the two sur- 
face fluxes may be required to produce the large vertical 
velocities associated with D. 

In summary, the surface fluxes of sensible heat and mois- 
ture contribute to a reduction in low-level static stability 
in the diabatic model atmosphere. This reduction is suffi- 
cient to shift the wavelength of maximum growth in the model 
atmosphere to shorter wavelengths than those found in the 
adiabatic model atmosphere. In addition, latent heating in 
the diabatic experiments enhances the development of shorter 
wavelengths. This latent heating is very sensitive to the 
availability of low-level moisture and, to a lesser degree, 
surface sensible heat fluxes. The fall case with a weaker 
jet and higher sea-surface temperatures is much more sensi- 
tive to changes in the surface fluxes. 
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V. SECONDARY LOW DEVELOPMENT EXPERIMENT 



Several observational studies of secondary lows or polar 
lows have been performed. Nitta and Yamamoto (1974) and 
others investigated the development of weak, small-scale 
cyclones over the East China Sea. Mullen (1979) and Reed 
(1979) have completed case studies of polar lows which form 
in polar outbreaks associated with major frontal systems or 
large-scale cyclones over the eastern Pacific Ocean. Rasmussen 
(1979) has investigated small-scale, intense cyclones in the 
North Atlantic Ocean which he also terms polar lows. These 
systems have two common features. They are small or inter- 
mediate scale (zonal wavelength of 1000 to 2000 km) , and they 
originate over the ocean. Reed (1979) has also noted that very 
few intermediate-scale cyclones occur over land. Because 
these lows appear mainly over the ocean, it is felt that air- 
sea fluxes are likely to be important during their evolution. 

Intermediate-scale waves are present throughout the dia- 
betic model experiments, although they are dominated by the 
large-scale primary cyclone. This chapter examines the evo- 
lution of the secondary cyclones using the technique described 
above of selectively removing the surface flux processes from 
the diabatic model. Table 5-1 is a list of the integrations 
performed for this experiment. Fields from the full diabatic 
experiment are used to initialize the other model experiments 
at day 4 (S4) and day 9 (S9) . 
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Table 5-1 



Experiments for the secondary development study 

* * 
Exp . Initial Integration Model 





Fields 


Period 


Version 


D 


S 


(wave no. 6) 


15 days 


complete diabatic 


D-FS(4) 


S4 

at 


(fields from D 
day 4) 


6 days 


diabatic w/o surface 
sensible heat flux 


D-FQ (4) 


S4 




6 days 


diabatic w/o surface 
moisture flux 


D-FSQ(4) 


S4 




6 days 


diabatic w/o surface 
moisture or sensible 
heat fluxes 


F(4) 


S4 




6 days 


adiabatic w/friction 


D-FS (9) 


S9 

at 


(fields from D 
day 9) 


3 days 


diabatic w/o surface 
sensible heat flux 


D-FQ (9) 


S9 




3 days 


diabatic w/o surface 
moisture flux 


D-FSQ(9) 


S9 




3 days 


diabatic w/o surface 
moisture or sensible 
heat fluxes 


F(9) 


S9 




3 days 


adiabatic w/friction 



See Sections II. C and II. D for a more complete dis- 
cussion of model variations and initial conditions. 
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Four cases involving development of a secondary cyclone 
wave are examined. Surface pressure maps for an illustra- 
tive time from each case are presented in Figs. 5-1 a-d. 

The first three maps are from D and the last is from D-FSQ. 

A common feature of the first three cases is that the second- 
ary lows develop in the equatorward flow over warmer water, 
similar to the polar lows studied by Reed (1979) and Rasmussen 
(1979) . The first two cases examined (see Figs. 5-la and b) 
are fall hemisphere developments, while the last two cases 
(Figs. 5-lc and d) are spring hemisphere developments. 
Intermediate-scale waves grow more readily in the initial de- 
velopment experiments. Not only is the number of secondary 
lows larger in the fall hemisphere than in the spring hemis- 
phere of the diabatic experiment, they are also more intense. 
It will be shown that the secondary lows in the fall case may 
evolve into the next primary cyclone. In the spring case, 
the secondary lows remain relatively weak compared to the 
primary cyclone and dissipate more rapidly. 

Expanded analyes of various model atmosphere parameters 
and derived quantities are presented below. Analyses cover 
an area of 18° lat by 22.5° long or 10 x 10 model grid points. 
The expansion is about the same center point for any given 
time and hemisphere, and is chosen as approximately the posi- 
tion of the center of the secondary cyclone in D-FS, except 
as noted. Southern hemisphere analyses are inverted with 
south at the top of the map for ease of comparison with the 
northern hemisphere cases. 
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Figure 5-1. Surface pressure maps for D for a. Case A, 
day 6.5, b. Case B, day 11, c. Case C, 
day 11 and d. Case D, D-FSQ, day 10.5. 
Contour interval 4 mb. Secondary low for 
case study indicated by an arrow in each 
figure . 
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A. CASE A- -FALL HEMISPHERE DAYS 6-8 



Case A gives an example of selective development of a 
secondary low. At day 6 of D a double low exists approximately 
35 degrees west of the main primary low (see position in Fig. 
3-2a) . A very similar system exists in the same location in 
experiments D-FS, D-FQ and D-FSQ (see Figs. 3-7 a-c, respec- 
tively) . At day 6 all four diabatic experiments, D, D-FS, 

D-FQ and D-FSQ, exhibit a dual low structure with a 1004 mb 
central pressure (see expanded surface pressure analyses 
Figs. 5-2 a-d) . These secondary lows evolve differently dur- 
ing the next 12 h (Figs. 5-3 a-d and 5-4 a-d). The cyclones 
in D-FS and D-FSQ (no surface sensible heat flux) are more 
intense than their counterparts, D and D-FQ. In addition, 
in D-FQ and D-FSQ (no surface moisture flux) the westward 
lobe of the initial dual low (see Fig. 5-2a) has developed, 
causing a phase lag of 9° long behind the lows in D and D-FS 
after 12 h. A slight phase lag (1° long) is also associated 
with the addition of the surface sensible heat flux. The low 
in D lags that in D-FS and, similarly, the low in D-FQ lags 
that in D-FSQ. 

In D and D-FS, the secondary cyclone intensifies as it 
moves into the northerly flow region of the primary cyclone 
(see Figs. 5-5 a-d) , and eventually absorbs and replaces the 
primary cyclone. This is not true in D-FQ and D-FSQ. In 
the latter experiments the secondary cyclone gradually weakens 
and is absorbed by the primary cyclone. The secondary lows 
are relatively sensitive to the air-sea fluxes. Therefore, 
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Figure 5-3. As in Fig. 5-2 except for day 6.25. 
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modification of the fluxes can result in a greatly different 
evolution in the atmosphere on a time scale of three to four 
days. In this case, the large differences are due more to 
selection of a different perturbation, rather than being due 
to significant differences in intensification rates. 

The major differences in development occur in response to 
the low-level temperature advection and the release of latent 
heat. The low develops within a strong baroclinic zone (see 
Figs. 5-6 a-d) with the strongest temperature gradient and 
thermal advection located to the east of the low. The low- 
level temperature gradient is largest in D-FSQ (no surface 
flux) and smallest in D (full diabatic) . 

Precipitation is restricted to the region east of the 
secondary lows in all of the experiments (see Figs. 5-7 a-b 
and 5-8 a-d). At day 6 (Fig. 5-7), the precipitation rate 
associated with the developing secondary low in D and D-FS is 
up to 11 cm/day, although this rate decreases slightly at day 
6.25 (Fig. 5-8). The precipiation rate in D-FQ and D-FSQ 
is less than 2 cm/day at day 6 (not shown) and increases 
slightly to 2 and 5 cm/day, respectively, at day 6.25. The 
maximum precipitation occurs over the region of maximum 
horizontal temperature advection (Fig. 5-6) in all experiments 
Convective precipitation is less than 20% of the total precipi 
tation in all cases. The low-level specific humidity fields 
at day 6 (Fig. 5-9 a-d) reflect the loss of the surface 
moisture flux in D-FSQ and D-FQ. Specific humidity is approxi 
mately 30% higher in the region of maximum precipitation in 
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Figure 5-7. Precipitation rate for day 6 and expanded 
region as in Fig. 5-2 for a. D and b. D 
Contour interval 2 cm/day. Experiments D 
and D-FSQ not shown as rate less than 2 
cm/day . 
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Figure 5-9. Specific humidity at the lowest model sigma level (900 
mb) for day 6 and expanded region as in Fig. 5-2 for 
a. D, b. D-FS, c. D-FQ and d. D-FSQ. Contour 
interval 1 g/kg . 




D and D-FS. The pattern of maximum specific humidity, and 
also the low-level wind field seems to indicate strong mois- 
ture advection from the southwest of the secondary low. In 
D-FQ and D-FSQ only, a region of relatively low specific 
humidity is present in the western quadrant of the secondary 
low. 

The direction of motion of the secondary low is toward 
the southeast in all of the experiments. This is toward the 
maximum low-level thermal gradient, the area of highest precipi- 
tation and largest low-level moisture. The implication is 
that the development is brought about by low-level thermal 
processes as noted by Petterssen a^. , (19 62), and that the 

surface fluxes of sensible heat and moisture significantly 
modify this development. 

To investigate the relative importance of the surface 
fluxes in the development of the secondary cyclone, the ther- 
mal terms of Petterssen 's development equation (Petterssen, 

1956) are evaluated. Petterssen 's development equation de- 
rived in a sigma coordinate system is: 



^^gO _ ^^gu _ g „2 8h 

at at “ f at ' 



A B C 



(5-1) 



where: 

C = geostrophic vorticity 

g = acceleration of gravity 
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f = Coriolis parameter, and 



h 



thickness from level 0 to level u. 



Here, term A represents the local change in geostrophic 
vorticity at the surface, term B is the local change in geo- 
strophic vorticity at an upper level and term C is the change 
in thermal vorticity due to a change in thickness caused by 
warming or cooling the column of air. 

Terms A and B can be inferred from the surface pressure 
fields and upper-level wind fields. Term C can be expanded 
to evaluate the contributions of the individual thermodynamic 
processes as follows: 





D 



E 




( 5 - 2 ) 



P 



F 



where : 



V 



H 



horizontal wind on a sigma surface 



V 



a 



horizontal gradient operator on a sigma 
surface 



T 



ir 



temperature 
surface pressure 



a = 


specific volvime 








C 

P 


specific heat at 


constant pressure 




a = 


P/tt where p is pressure on a 


sigma surface 


a = 


da/dt 








dH/dt = 


heating rate due 
as radiation and 


to diabatic 
latent heat 


effects 
release , 


such 

and 


R 


gas constant for 


dry air. 







Here, term D represents the change in thickness due to tempera- 
ture advection, term E the change due to dynamic warming or 
cooling and term F the change due to diabatic heating. In 
general, warm advection and precipitation (terms D and F) 
tend to enhance low-level vorticity, whereas the compensating 
upward vertical motion (term E) tends to oppose development. 

Certain assumptions are necessary to evaluate the terms 
of the above equations. The upper level is chosen as sigma 
level 4 which is approximately 500 mb. Only a very weak dis- 
turbance in essentially straight flow exists at this level. 

The equations were also evaluated at one of the time periods 
using sigma level 5 as the upper level (approximately 700 mb) 
and yielded qualitatively the same results. The lowest model 
layer contains an idealized planetary boundary layer. Cal- 
culations in this layer are performed using an average of 
the sigma level 6 (900 mb) values and the surface values (19.5 
m for wind values) , rather than recalculate the complete 
boundary layer formulation . 
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Term B of the Petterssen development equation represents 
the contribution to low-level vorticity due to the change in 
upper-level vorticity. Comparison of sigma level 4 (500 mb) 
winds from D and D-FQ at days 5.5, 6 and 6.25 revealed nearly 
identical (± 1 m/s) winds with a small dual trough located 
vertically above the surface secondary lows in both experi- 
ments. The contribution to the secondary development from 
term B appears to be minor and in any case should be similar 
in the two experiments. 

Term D represents the change in thickness due to horizon- 
tal temperature advection. The region of strongest advection 
is discussed above in terms of the low-level temperature and 
wind fields (Fig. 5-6) . Term D (see Figs. 5-10 a-d) gives 
the thickness change in the lower model atmosphere produced 
by the horizontal temperature advection. The smallest thick- 
ness increases east of the secondary lows occur in D, and 
the largest in D-FSQ. The area of positive thickness increases 
is also much larger in D-FS, D-FQ and D-FSQ (no surface flux) . 
This maximum thickness increase occurs five degrees westward 
in D-FQ and D-FSQ (no surface moisture flux experiments) . 

This is consistent with the development of the western per- 
turbation in these experiments. 

Term E (Figs. 5-11 a-d) represents the dynamic cooling 
of the rising air and acts as a braking term to oppose develop- 
ment. This term is influenced both by the lapse rate of the 
environmental air and by the magnitude of the vertical motion — 
either larger upward motion or a smaller lapse rate leads to 
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Figure 5-11. As in Fig. 5-10 except for term 



greater braking. D and D-FS have very large values of term 
E. The secondary low in D has the highest central pressure 
at day 6.25 (see Fig. 5-3) , and yet the magnitude of term E 
is twice that of D-FQ and D-FSQ. As discussed in the initial 
growth experiments, this is associated with vertical veloci- 
ties that are two to three times larger in D and D-FS . 

Term F (Fig. 5-12 a-d) represents the warming/cooling of 
the column of air through the diabatic processes of radiation, 
precipitation (latent heat release) and surface sensible heat 
flux. The greater warming in D-FS relative to D seems to be 
related to the 20% greater precipitation (see Fig. 5-7) in 
D-FS at day 6.25. This difference obscures the influence of 
the surface sensible heat flux between D and D-FS. 

Diabatic heating (term F) contributes more to the develop- 
ment of the secondary low than does thermal advection (term 
D) . However the mechanism for lifting the air and causing 
the diabatic heating is clearly the thermal advection. This 
is because of the near colocation of the maxima in these fields 
in Fig. 5-10 and Fig. 5-12 (also see Figs. 5-6 and 5-8). The 
loss of the surface fluxes of sensible heat and moisture lead 
to a larger contribution from thermal advection. The loss of 
the surface moisture flux, however, has a much greater effect, 
because it influences both the location of phasing as well as 
the intensity of the development. The total contribution of 
the thermal processes (term C) to the thermal vorticity (Fig. 
5-13) is several times larger in the experiments with surface 
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Figure 5-12. As in Fig. 5-10 except for term 
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Figure 5-13. As in Fig. 5-10 except for term C. Contour interval 



moisture flux. This process appears to be the reason why 
the secondary lows continue to develop in D and D-FS while 
they slowly dissipate in D-FQ and D-FSQ. 

Although there is some deep convection, precipitation 
caused by large-scale lifting appears to be responsible for 
the development of the secondary lows in D and D-FS. This 
lifting is caused by warm advection and occurs in a region 
of strong moisture gradient. This moisture gradient is main- 
tained by the advection of moisture and by the surface mois- 
ture flux and results in much greater precipitation in D and 
D-FS. 

B. CASE B — FALL HEMISPHERE DAYS 10-12 

The secondary development in this case is similar to the 
one examined in Case A. Model versions D-FS, D-FQ and D-FSQ 
are also initialized with fields from D at day 9. A 1008 mb 
secondary low first appears at approximately day 10 in all 
model versions. At day 10.5, 1.5 days after initialization, 
the secondary low enters the southwest quadrant of the primary 
low (see Fig. 5-14 a-d and also Fig. 5-lb) with approximately 
the same central pressure, 1007 mb ± 1 mb, and location in 
all the experiments. During the next 18 h, the secondary low 
develops (see Figs. 5-15 a-d and 5-16 a-d) much more rapidly 
in D and D-FS, and becomes the primary low by day 12. In 
D-FQ and D-FSQ (no surface moisture flux) the secondary low 
slowly dissipates after day 11. Also the phase speed of the 
secondary low is less in D-FQ and D-FSQ. This results in a 



103 






^ CN 

^ Q 
O -H 

r— I • fd 

> 

>1 M 
U 0) 
(TJ M -P 
O C 

C M-t -H 

O 

U O P 
<D P 13 
U) (1) O 
^ -P 
<D C 
^ U) O 
-P -H U 

s 

m 0) 
o x: • 
a 

U) iH CO 

cij «H (J 4 

<d I 
g Q 

CD 0) • 

p x: 

-p 

0 ) 

CO C T3 
0) -H C 
P (TJ 

am 

• O 

0 ) o a 

U rH I 
03 Q 
^ >1 
P rrJ • 
P T5 U 
CO 

T5 fO ^ 
0) CO 
T5 CQ a 
C I 
<T3 0) Q 
a CO 

w CJ XJ 



r-l 

I 

m 



0 ) 

p 

:3 

CP 

•H 

a 



104 




105 




106 



Figure 5-16. As in Fig. 5-14 except for day 11.25. 



phase lag of 6° long relative to the secondary low in D during 
the period from day 10.5 to day 11.25. 

Case B closely parallels Case A. To simplify the dis- 
cussion, only the contributions of terms D, E, and F in the 
Petterssen development equation will be presented. Term D 
(Fig. 5-17 a-d) yields similar contributions as in Case A. 

D-FS and D-FSQ (no surface sensible heat flux) , have larger 
contributions to development than the experiments with sensible 
heat flux, D and D-FQ, respectively. Term E (Fig. 5-18 a-d) 
again represents a much larger braking in D and D-FS due to 
the larger vertical motions associated with the secondary low. 
The diabatic heating term (term F, Fig. 5-19 a-d) , represents 
the largest contribution to the thickness increase. The sum 
of terms D and F is not offset by term E in D and D-FS, which 
allows the secondary low to develop in D and D-FS, but dissi- 
pate in D-FQ and D-FSQ. The thermal advection term contributes 
to larger development in D-FS relative to D and less dissi- 
pation in D-FSQ than in D-FQ. Lifting due to warm advection 
is clearly the dominant mechanism in producing precipitation, 
and thermal advection is promoted by the absence of the sur- 
face sensible heat flux. 

C. CASE C — SPRING HEMISPHERE DAYS 10-12 

Case C is an example of a secondary low that fails to 
deepen. As in Case B, model versions D-FS, D-FQ and D-FSQ 
are initialized with fields from D at day 9. A secondary low 
is evident in all model versions about day 10.5, or 1.5 days 
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C LONGITUDE ^ LONGITUDE 

Figure 5-17. Term D of the Petterssen development equation for day 11 

and an expanded region as in Fig. 5-15 for a. D. b. D-FS 
c. D-FQ and d. D-FSQ. Contour interval 5 x 10“^ °C/s. 
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after initialization. The low is in approximately the same 
location (see Fig. 5-20 a-d and Fig. 5-lc) with only small 
differences in initial intensity. The secondary low moves 
through the region southwest (would be northeast if plotted 
normally) of the primary low, which corresponds to the loca- 
tions of the developments that occurred in Cases A and B. In 
this case, however, the low merely maintains the same inten- 
sity for approximately 12 h, and then slowly dissipates in 
all model versions (Figs. 5-21 and 5-22) . 

The Petterssen development equation terms have been evalu- 
ated for Case C at day 10.75 (Figs. 5-23, 5-24 and 5-25). 

Term D (Fig. 5-23) reflects a somewhat weaker thermal advec- 
tion contribution than in Cases A and B, although the decrease 
in Case C is less than 10% for D and D-FS . It is the second- 
ary lows in these two model versions which develop strongly 
in Cases A and B. Term F (Fig. 5-25) demonstrates a decrease 
of 60% to 80% in the diabatic warming of Case C from Cases A 
and B (Figs. 5-12 and 5-19) in experiments D and D-FS. The 
maximum values of the diabatic heating in Case C are only of 
the same order as those in D-FQ and D-FSQ of Cases A and B. 

The secondary low in D-FQ and D-FSQ dissipated eventually in 
both Cases A and B. The implication here is that a large 
contribution from diabatic heating, term F, is required for 
the secondary low to continue to develop . 

D. CASE D — SPRING HEMISPHERE DAYS 9-11 

This case differs from the previous three cases in that 
the secondary low develops in the poleward flow in the eastern 
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Figure 5-21. As in Fig. 5-20 except for day 10.75. 
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Figure 5-23. Term D of the Petterssen development equation for day 10.75 
and an expanded region as in Fig. 5-21 for a. D, b. D-FS, 
c. D-FQ and d. D-FSQ. Contour interval 5 x 10“5 °c/s. 
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quadrant of the primary cyclone (see Fig. 5-ld) . Also, in 
apparent contradiction to the first three cases and the re- 
sults discussed in Chapter III, the largest development occurs 
in the adiabatic experiment, F, with the least development 
in D. 

Figs. 5-26 a-b present the expanded surface pressure maps 
for F corresponding to those presented in Figs. 5-20 a-d and 
5-22 a-d for the diabatic experiments. The secondary low 
reaches its maximum intensity between days 10.5 and 11 in 
all experiments, prior to entering the region of strong 
equatorward flow where development occurs in Cases A and B. 

At this time, day 10.5, the secondary low is deepest in F, 
followed by D-FSQ, and with the highest central pressure 
occurring in D. 

Temperature fields for experiments D and F (see Figs. 

5-27 a-b) reveal a much larger low-level temperature gradi- 
ent to the east of the primary cyclone in F than D. The 
gradients for D-FSQ, D-FS and D-FQ (not shown), respectively, 
are also larger than D and the secondary lows are corres- 
pondingly more intense (see Figs. 5-20 a-d). Expanded plots 
of the low-level temperature fields (Figs. 5-28 a-b) reveal 
a temperature gradient in F approximately twice that in D. 
Thus, the indication is that low-level temperature advection 
is responsible for the difference in development in this case. 

In terms of the Petterssen development equation, the 
secondary low developing in F (adiabatic) has no contribu- 
tion from diabatic heating (term F) , yet it develops more 
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Figure 5-26. Expanded surface pressure maps of secondary 

low in F in the spring hemisphere for a. day 
10.5 and b. day 11. Region as in Figs. 5-20 
and 5-22, respectively. Contour interval 
2 mb. 




Figure 5-27. Southern hemisphere temperature fields at 
sigma level 6 (900 mb) at day 10.5 for 
a. D and b. F. Contour interval 5°C. 
Location of secondary low indicated by X. 
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Figure 5-28. 



As in Fig. 5-27 except for expanded region 
as in Fig. 5-26a. Contour interval 2°C. 




Figure 5-29. Term D of the Petterssen development equa- 
tion for day 10.5 and an expanded region 
as in Fig. 5-26a for a. D and b. F. 
Contour interval 5 x 10“5 «c/s. 
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than the secondary low in D. Figs. 5-29 a-b give the con- 
tribution of the low-level temperature advection (term D) 
to development of the secondary low at day 10.5. The region 
of positive thickness change (negative values in figure) is 
larger as well as more intense in F. 

Developments of this type (Case D) in the eastern quadrant 
of the primary cyclone occur in both hemispheres, although 
the intensity changes are generally not as marked as in this 
example. These situations, in which the air-sea fluxes and 
diabatic effects appear to hinder the cyclone development, 
are more common and more marked in the coarse-mesh model 
experiments. It appears that the relative role of the physi- 
cal processes is changed in the different resolutions. In 
the northern or fall hemisphere, the larger diabatic effects 
in experiments D and D-FS tend to mask the contribution of 
the stronger low-level temperature advection. 

E . SUMMARY 

These cases clearly demonstrate the sensitivity of the 
model development of secondary lows to both the surface fluxes 
of sensible heat and moisture and to the baroclinicity of the 
model atmosphere. The loss of the surface sensible heat flux 
contributes to a greater low-level temperature gradient 
in both the spring and fall hemispheres. This contributes to 
an increase in low-level temperature advection, and leads to 
a greater development in the experiments with no sensible heat 
flux. The surface moisture flux is responsible for much 
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larger differences in the development of secondary lows. 

Here the physical mechanism is the latent heat release 
triggered by vertical velocities induced by the low-level 
temperature advection. 

In the fall hemisphere with the weaker jet, the secondary 
low develops rapidly in the experiments with a surface mois- 
ture flux, and eventually replaces the primary, large-scale 
cyclone. In the experiments without surface moisture flux, 
the secondary lows slowly dissipate. In addition, there is 
a six-to-ten degree phase lag of the secondary lows in these 
experiments relative to those which had surface moisture 
fluxes. This phase lag occurs primarily during the 12 to 
18 h period in which the secondary lows develop most rapidly. 

In sharp contrast, in the spring hemisphere with the 
stronger jet and lower sea-surface temperatures, the second- 
ary lows dissipated after one or two days, without intensi- 
fying greatly. The low-level temperature advection in this 
hemisphere is only about 10% less than in the experiments in 
which the secondary lows develop strongly in the fall hemis- 
phere. The contribution from the diabatic heating, however, 
is much less in the spring hemisphere experiments with sur- 
face moisture flux, and is comparable to the experiments in 
which the secondary low dissipated in Cases A and B. 

Reed (1979) concluded that the polar lows in the eastern 
North Pacific Ocean were predominately a baroclinic phenomena 
with possible contributions from thermal and barotropic proc- 
esses. Rasmussen (1979), in a similar study of eastern North 
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Atlantic Ocean polar lows, concluded that they were pre- 
dominately a thermal instability phenomena. The results of 
the model experiments presented here demonstrate that both of 
the above situations could exist in the marine atmosphere. 

The instability process which dominates in any given instance 
is sensitive to the surface fluxes of sensible heat and mois- 
ture, as well as to the degree of baroclinic instability 
present in the atmosphere. In these experiments, the surface 
fluxes are critically important on a time scale of two to 
four days, which is well within the predictability of the 
real atmosphere. 
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VI. CONCLUSIONS 



The objective of this research is to use a numerical model 
to examine the role of the air-sea fluxes in extratropical 
cyclogenesis over the open ocean. The method employed is to 
specify analytic initial atmosphere and ocean conditions which 
are typical of fall and spring conditions in the real atmos- 
phere. A state-of-the-art atmospheric prediction model is 
used to examine the contributions of the surface fluxes of 
moisture and sensible heat. This is accomplished by selec- 
tively removing the surface flux parameterizations from the 
model and comparing the results with a control run generated 
with the complete diabatic model. Consequently, this re- 
search emphasizes the role of the individual processes in the 
model and in the atmosphere rather than the response of the 
model to small changes in the sea-surface temperatures (e.g., 
Spar and Atlas, 1975) . 

The surface fluxes make a significant contribution in 
these simulations of cyclone evolution. In the fall hemi- 
sphere with a weaker upper-level jet and higher sea-surface 
temperature, the removal of the surface fluxes results in a 
large reduction in the initial growth rate of the model cy- 
clones. In the presence of the more intense jet and lower 
sea-surface temperature of the spring hemisphere, removal of 
the surface fluxes has a smaller impact on the model cyclone 
growth. The contribution of the surface fluxes to the growth 
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of a pre-existing, large-scale cyclone is found to be rela- 
tively small as expected from previous studies. 

The smaller-scale secondary cyclones respond rapidly to 
the removal of the model surface fluxes. Removal of the sur- 
face moisture flux prevents the development of secondary 
lows which might otherwise eventually replace the primary 
cyclone, as occurs in the complete diabatic model. Removal 
of the surface sensible heat flux results in the development 
of more intense cyclones although the overall evolution of 
the cyclone is not altered from that in the complete dia- 
batic experiment. The result that different perturbations 
eventually become the primary cyclone after three to four days 
is significant. One may conclude that extending the useful- 
ness of numerical model forecasts of open-ocean cyclones 
beyond the present limit of two to three days will require 
accurate specification of the air-sea fluxes. 

The experimental technique of removing the surface-flux 
parameterizations from the model has no real world analog. 

The modelled results can not be verified with observations 
and indicate only the sensitivity of the model cyclones to the 
air-sea fluxes. The ability to perform experiments within a 
controlled, but admittedly idealized, environment can provide 
insights, however, that would otherwise require very large, 
detailed statistical studies. These results do suggest why 
some disturbances over the ocean develop while other, seem- 
ingly similar, disturbances do not develop. 
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The surface fluxes of moisture and sensible heat act to 



reduce the large gradients of low-level temperature and mois- 
ture in the model atmosphere. These gradients occur in re- 
gions of strong equatorward (poleward) advection associated 
with the primary cyclone, especially over the zone of largest 
sea-surface temperature gradient. Large-scale lifting produced 
by the low-level temperature advection triggers the release 
of latent heat as hypothesized by Petterssen et , (1962). 

These two processes jointly lead to a large increase in ver- 
tical motion and, therefore, cyclone development. Removal of 
the surface sensible heat flux allows a larger low-level 
temperature gradient, resulting in enhanced vertical motions 
and in increased development of the cyclones. Removal of 
the surface moisture flux results in a large decrease in pre- 
cipitation or latent heat release and, therefore, a large 
decrease in cyclone development. 

Previous studies have focused on the region behind, or 
westward, of the cyclone, where the strongest fluxes occur. 

The cooperative effect of lifting produced by low-level tempera- 
ture advection and the release of latent heat, however, results 
in the largest cyclone response occurring ahead, or eastward, 
of the cyclone. The results of these experiments indicate 
that the magnitude of the surface fluxes in terms of the 
amount of energy transferred across the air-sea interface is 
not as important as the location in which they occur. This 
resolves the apparent contradiction found by Wei (1979) and 
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also implies that accurate specification of the sea-surface 
temperature distribution is crucial to the evolution of open- 
ocean cyclones, at least on a medium range time scale of 
three to five days. 

This experiment addresses open-ocean conditions with 
fixed sea-surface temperatures and moderate values of surface 
moisture and sensible heat fluxes. A variety of cyclone evo- 
lutions similar to those observed by Nitta and Yamamoto (1974) , 
Reed (1979) and others is produced in the model integrations. 
Much more intense over-ocean cyclogenesis has been observed 
(Sanders and Gyakum, 1980; Pyke, 1965) in association with 
land-sea and ice-sea boundaries where the atmosphere is far 
from equilibrium with the ocean surface fluxes. Further re- 
search is necessary to investigate more fully the role of 
the surface fluxes in these cases of more intense cyclogene- 
sis. Numerical experiments of the type described here can 
be helpful in understanding the physical processes involved 
in these cases, and will hopefully lead to improved prediction 
of extratropical cyclones over the oceans . 
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